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The optical excitations in C70 and higher fullerenes, including isomers of 

>< 

' C76, C78, and C84, are theoretically investigated. We use a tight binding 



model with long-range Coulomb interactions, treated by the Hartree-Fock and 
configuration-interaction methods. We find that the optical excitations in the 
energy region smaller than about 4 eV have most of their amplitudes at the 
pentagons. The oscillator strengths of projected absorption almost accord with 
those of the total absorption. When the projection is performed on each pen- 
tagon and pentagon dimers, the resultant spectrum in the low energy region 
is quite different from that of the total absorption. The spectral shapes of the 
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total absorption are turned out to be determined mainly by the geometrical 
distributions of the pentagons in the fuUerene structures. 
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- 2 - 



1. Introduction 



Recently, the fuUerene family Cn with hollow cage structures has been inten- 
sively investigated. A lot of optical experiments have been performed, and 
excitation properties due to 7r-electrons delocalized on molecular surfaces have 
been measured. For example, the optical absorption spectra of Ceo and C70 
[1,2] have been reported, and the large optical nonhnearity of Cgo [3,4] has 
been found. The absorption spectra of higher fuUerenes (Cyg, C78, C84, etc.) 
have also been obtained [5,6]. For theoretical studies, we have apphed a tight 
binding model [7] to Cgo, and have analyzed the nonlinear optical properties. 
Coulomb interaction effects on the absorption spectra and the optical nonhn- 
earity have been also studied [8]. We have found that the hnear absorption 
spectra of Ceo and C70 are well explained by the Frenkel exciton picture [9] 
except for the charge transfer exciton feature around the excitation energy 2.8 
eV of the Ceo solids [2] . Coulomb interaction effects reduce the magnitude of 
the optical nonhnearity from that of the free electron calculation [8] , and thus 
the intermolecular interaction effects have turned out to be important. 

In the previous paper [10], we have extended the calculation of Ceo [9] to 
one of the higher fuUerenes Cye. We have discussed variations of the opti- 
cal spectral shape in relation to the symmetry reduction from Ceo and C70 
to Cre' the optical gap decreases and the spectra exhibit a larger number of 
small structures in the dependences on the excitation energy. These prop- 
erties seem to be natural when we take into account of the complex surface 
patterns composed of pentagons and hexagons. In order to understand the 
patterns clearly, the idea of the "phason hne" has been introduced [11] using 
the projection method on the honeycomb lattice plane [12]. There are twelve 
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pentagons in Cye- Six of them cluster on the honeycomb lattice, with one 
hexagon between the neighboring two pentagons. There are two groups of the 
clustered pentagons. The "phason hne" runs as if it divides the two groups. 

The purpose of this paper is to investigate relations between the optical 
properties and geometric structures in higher fuUerenes. The phason lines can 
characterize the geometries of fuUerenes. We define the "pentagonal" carbons 
as the atoms located at the points of pentagons. Then, the carbon atoms along 
the phason lines can be regarded as the remaining atoms after the petangonal 
carbons are taken out. We will calculate the optical absorption spectra at a 
certain combination of the pentagonal carbons. We project wavefunctions to 
selected pentagons, and calculate dipole moments using the projected wave- 
functions. Thus, the contributions from a part of the fuUerene to the optical 
spectra can be extracted. 

The main results are as follows: 

(1) The optical excitations in the energy region smaller than about 4 eV have 
most of their amplitudes at the pentagonal carbons. The oscillator strengths 
of absorption projected onto these carbons almost accord with those of the 
total absorption. 

(2) When the projection is performed on a smaller set of the pentagonal car- 
bons, the resultant spectrum in the low energy region is quite different from 
that of the total absorption. The structures in the absorption cannot be de- 
composed into contributions from subsets of pentagonal carbons. 

(3) The contributions from pentagon dimers in several kinds of higher fuUerenes 
(e.g., two isomers of Cyg and C84) are compared. We find that small peak struc- 
tures are mutually different. The spectral shapes of the total absorption are 
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mainly determined by the geometrical distributions of the pentagons in the 
fuUerene structures. 

In the next section, our tight-binding model is introduced and the calcu- 
lation method is explained. Sections 3 and 4 are devoted to numerical results 
and discussion about the properties of the absorption spectra. The paper is 
closed with a summary in section 5. 



2. Model 

We use the following Hamiltonian: 

H^Ho + H,^t. (1) 
The first term of eq. (1) is the tight binding model: 

Ho^-tYl i4,aCj,a + h.C.), (2) 

where t is the hopping integral and q^o- is an annihilation operator of a tt- 
electron with spin a at the ith carbon atom of the fuUerene. It is assumed 
that t does not depend on the bond length, because main contributions come 
from excitonic effects due to the strong Coulomb potential. The results do 
not change so largely if we consider changes of hopping integrals by bond 
distortions. We assume the following form of Coulomb interactions among 
TT-electrons: 

+ E^nn,,)(Ec.W-i)(EcW-i), (3) 

i^j a- T 

where r^j- is the distance between the ith and jth atoms and 

W{r) = ^ (4) 

^{l/Uf + {r/roVf 



is the Ohno potential. The quantity U is the strength of the on-site interaction, 
V means the strength of the long-range Coulomb interaction, and Tq is the 
average bond length. 

The model is treated by the Hartree-Fock approximation and the single 
excitation configuration interaction method, as was used in the previous pa- 
pers [9,10]. In ref. 9, we have varied the parameters of the Coulomb inter- 
actions, and have searched for the data which reproduce overall features of 
experiments of Ceo and C70 in solutions. We have found that the common 
parameters, U = it and V — 2t, are reasonable. Thus, we use the same pa- 
rameter set for higher fuUerenes. The quantity t is about 2 eV as shown in ref. 
9. The Coulomb interaction strengths depend on the carbon positions. We use 
the lattice coordinates which are obtained by the pubhc program "FULLER" 
[13,14]. The optical spectra become anisotropic with respect to the orienta- 
tion of the molecule against the electric field of light, as reported in the free 
electron model (Hiickel theory) [15]. We obtain numerical optical absorption 
spectra by summing the data of three cases, where the electric field of light is 
along the x-, y-, and z-axes. 

We use a projection operator in order to extract contributions to the optical 
spectra from a certain part of fuUerenes. If we write the a projection operator 
to a part of lattice site set as P, the oscillator strength between the ground 
state \g) and the excited state \k) is written: 




E,[\{K\PxP\g)\'' 



+ 



- P)x{l - P)\g)\' 



+ 



{g\PxP\K){K\{l-P)xil-P)\g) 
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+ {g\{l-P)x{l-P)\K){K\PxP\g)l 



(5) 



where is the excitation energy, and the electric field is parallel with the 
X-axis. In eq. (5), the first term is the contribution from the projected part, 
and the three other terms are the remaining part. The total optical absorption 
is calculated by the formula: 



where p{ijj) — ^ /[k{uP -\-^'^)\ is the Lorentzian distribution of the width 7. The 
projected absorption is calculated by eqs. (5) and (6). The projected part does 
not satisfy a sum rule. So, this results in a singularity where excitation energy 
is large. We will discuss the optical spectra in the energy region far from the 
singularity. 

3. Optical absorption in C70 and C76 

Figure 1 shows the molecular structures and optical spectra of the C70 molecule 
and C76 with the D2 symmetry, which have been found in experiments. The 
black atoms are the carbons along the phason lines. The hatched circles are 
the pentagonal carbons. In C70, the phason line runs along the ten carbons 
which are arrayed like a belt around the molecule. In Cye, the phason line is 
located almost along the outer edge of the molecule of Fig. 1(b). The total 
optical absorption is shown by the bold line, and the absorption from all the 
pentagonal carbons is shown by the thin line. We find that the optical excita- 
tions in the energy region lower than 2t are almost composed of the excitations 
at the pentagonal sites. This property is common to C70 and D2-C-jq, and also 
to the Td-CjQ for which the calculated data are not shown. In higher energy 
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regions, the thin hnes give relatively larger oscillator strengths, but this is 
an artifact of the projected wavefunctions. The absorption spectra calculated 
from the projected wavefunctions do not satisfy the sum rule, i.e., the area 
between the abscissa and the curve does not become constant regardless of 
the excitation wavefunctions. The similar artifact will be found in the figures 
shown afterwards. We believe that the projected optical absorption spectra 
are reliable in low energy regions only. Therefore, we limit our comparison of 
the spectra to the energy region lower than about 2i ~ 4 eV. 

In Cgo, the edges of the pentagons are the long bonds, and the bonds 
between the neighboring hexagons are short bonds. The wavefunctions of 
the fivefold degenerate highest-occupied-molecular-orbital (HOMO) have the 
bonding property, and that the threefold degenerate lowest-unoccupied-molecular- 
orbital (LUMO) has the antibonding property. As the carbon number in- 
creases, hexagons are inserted among pentagons. The wavefunctions near the 
LUMO of the higher fuUerenes still have the antibonding properties, thus they 
tend to have large amplitudes along the edges of pentagons which have the 
characters hke long bonds of Ceo- Recently, the bunching of the six energy lev- 
els higher than the LUMO has been discussed in the extracted higher fuUerenes 
[16] . The wavefunctions near the LUMO distribute on the pentagons. This fact 
can be understood as the properties characteristic to antibonding orbitals. As 
the excited electron mainly distributes at the pentagonal carbons, the electron- 
hole excitation has large amplitudes at these pentagons. Thus, the oscillator 
strengths of the low energy excitations are mainly determined by wavefunctions 
at the pentagonal carbons. This is the reason why the projected absorptions 
nearly accord with the total absorptions in the energy regions smaller than 
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about 2t. 

If the projections are performed onto each pentagon, we can know contri- 
butions to optical spectra from the projected carbon sites. We would like to 
look at this feature, for example, in D2-C7Q. There are three carbon atoms, 
which are not equivalent with respect to symmetries, in this isomer. These 
pentagons are indicated by the symbols, A-C, in Fig. 2(a). The projected 
absorption spectra are shown by thin curves, superposed with the total ab- 
sorption in Figs. 2 (b-d). The projected absorption is multiplied by the factor 
12, in order to compare with the total absorption. We find that the projected 
absorption exhibits small structures in the energy region smaller than 2t. The 
structures depend on the kind of carbons. The spectral shapes and oscillator 
strengths are much far from those of the total absorption. It would be diffi- 
cult to assign experimental features of the total absorption with a set of the 
limited number of carbon atoms. The excitation wavefunctions at the twelve 
pentagons give rise to the shape of the absorption spectra totally. 



4. Optical absorption in C78 and 

Figures 3 and 4 show the molecular structures and the calculated optical ab- 
sorption spectra of the extracted isomers of Cyg and C84, respectively. The 
notations in Figs. 3 and 4 are the same as in Fig. 1. The property, that 
the low energy excitations have most of their amplitudes at the twelve pen- 
tagons, is seen in these five isomers. The bold line and the thin line almost 
agree in the energy region lower than about 2i ~ 4 eV. It is of course that 
the agreement becomes a little bit worse from that in C70 and Cye of Fig. 1, 
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because the number of atoms along the phason hnes increases as the fuUerenes 
become higher. But, the similar wavefunction properties seem to persist in 
the calculated fuUerenes, C78 and C84. 

The pentagon dimers exist in the regions surrounded by the phason lines in 
the two C2V isomers of C78 [Figs. 3 (a) and (b)] and also in the isomer of C84 
[Fig. 4 (a)]. These three molecules have pentagon dimer structures commonly. 
The dimers are shown by the black pentagons in molecular structures of Fig. 
5. It is of some interests to look at whether these pentagon give rise to similar 
contributions to the optical spectra or not. The data of the thin lines are 
multiplied by the factor 6, in order to compare with the total absorption. We 
find that small peak structures are mutually different: for a typical example, 
the large feature around the energy l.bt represented by the thin fine in Fig. 5(c) 
is not present in Figs. 5(a) and (b). We thus conclude that the spectral shapes 
of the total absorption are mainly determined by the geometrical distributions 
of the pentagons in the fuUerene structures. 

5. Summary 

We have studied optical excitations in C70 and higher fuUerenes, including 
isomers of C76, C78, and C84. We have used a tight binding model with long- 
range Coulomb interactions, and have treated by the Hartree-Fock and CI 
methods. We have found that the optical excitations in the energy region 
smaller than about 4 eV have most of their amplitudes at the pentagonal 
carbons. Thus, the oscillator strengths of absorption projected onto these 
carbons almost accord with those of the total absorption. When the projection 
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is performed on a smaller set of the pentagonal carbons, for example, on each 
pentagon, the resultant spectrum in the low energy region is quite different 
from that of the total absorption. Therefore, the structures in the absorption 
cannot be decomposed into contributions from subsets of pentagonal carbons. 
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Figure Captions 



Fig. 1. Molecular structures and theoretical optical spectra of (a) D^d-Cro 
and (b) D2-C7Q. In the molecules, the black atoms are along the phason lines, 
and the hatched atoms are the pentagonal carbons. In the absorption spectra, 
the bold line is the total absorption, and the thin line is the absorption by the 
wavefunctions projected on the twelve pentagons. The units of the abscissa are 
taken as arbitrary, and the energy is scaled by t. The parameters are U — 4t, 
V = 2t, and 7 = OMt. 

Fig. 2. (a) The molecular structure of D2-C7G. The symbols, A-C, indicate the 
symmetry nonequivalent pentagons. The figures, (a), (b), and (c), compare 
the absorption projected on one of the three pentagons with the total absorp- 
tion. The bold line is the total absorption, and the thin line is the projected 
absorption. The units of the abscissa are taken as arbitrary, and the energy 
is scaled by t. The data of the thin line are multiplied by the factor 12. The 
parameters are U — At, V — 2t, and 7 = 0.06t. 

Fig. 3. Molecular structures and theoretical optical spectra of (a) C2V-C78, 
(b) the other C2^-C78, and (c) In the molecules, the black atoms are 

along the phason lines, and the hatched atoms are the pentagonal carbons. 
In the absorption spectra, the bold line is the total absorption, and the thin 
line is the absorption by the wavefunctions projected on the twelve pentagons. 
The units of the abscissa are taken as arbitrary, and the energy is scaled by t. 
The parameters are U = At, V = 2t, and 7 = 0.06t. 
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Fig. 4. Molecular structures and theoretical optical spectra of (a) D2d-Cs4 
and (b) D2-C84- In the molecules, the black atoms are along the phason lines, 
and the hatched atoms are the pentagonal carbons. In the absorption spectra, 
the bold line is the total absorption, and the thin line is the absorption by the 
wavefunctions projected on the twelve pentagons. The units of the abscissa are 
taken as arbitrary, and the energy is scaled by t. The parameters are U — 4t, 

V = 2t, and 7 = OMt. 

Fig. 5. Molecular structures and theoretical optical spectra of (a) C2^-C78, 
(b) the other C2v-C-j8, and (c) D2d-Cs4. In the molecules, the pentagon dimers 
are shown by the black pentagons. The each figure compares the absorption 
projected on the pentagon dimer with the total absorption. The bold line is 
the total absorption, and the thin line is the projected absorption. The units 
of the abscissa are taken as arbitrary, and the energy is scaled by t. The data 
of the thin line are multiplied by the factor 6. The parameters are U — 4t, 

V = 2t, and 7 = OMt. 
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